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Recent results on jets and the hadronic final state from the HERA collaborations HI and ZEUS are reviewed. 



1. INTRODUCTION 

Until the year 2000, the HERA experiments HI 
and ZEUS have collected integrated luminosities 
of about 130 pb -1 . A large number of interesting 
results on jets and the hadronic final state from 
this 'HERA 1' data taking period have already 
been published, and many more analyses are still 
ongoing. 

In 2001/2002 the HERA machine and also the 
experiments have undergone major changes with 
the aim of increasing the delivered luminosity by a 
factor of five. Due to technical problems, efficient 
data taking could not start before the end of 2003. 
Therefore, new results from the 'HERA 2' data 
taking period are only now starting to come out. 

In the recent years the analysis work on the 
HERA 1 data has concentrated on making op- 
timal use of the high statistics and either de- 
creasing statistical (and also systematic) errors or 
performing more and more differential measure- 
ments. In addition, increased statistics allows for 
the discovery or first measurement of new or ex- 
otic phenomena. 

This contribution will highlight new results 
concerning jets and the hadronic final state at 
HERA which were published or made preliminary 
since the EPS03 conference in Aachen. 

2. EVENT SHAPE VARIABLES 

Event shape variables like the thrust T, broad- 
ening B or the C parameter measure aspects of 
the topology of an event's hadronic hnal state. 
These variables are inclusive and infrared safe so 
that hadronisation corrections can be estimated 
using power correction (PC) models. 
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Figure 1. ZEUS event shape variables: Simulta- 
neous fit of as and the PC parameter oo to the 
mean of event shape variables using NLO+PC 
predictions. 



ZEUS has performed measurements of event 
shape variables in the Breit frame in the range 
80 < Q 2 < 20480 GeV 2 [T]. The choice of the 
Breit frame assures that the current region, i.e. 
the direction of flight of the scattered parton, re- 
sembles a single hemisphere in e + e~ collisions, 
thus making the HERA results comparable to re- 
sults from other colliders, especially LEP. 

The PC model suggests that non-perturbative 
corrections to the predictions of observables can 
be parametrised with a universal parameter 5q 
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according to (F) = (F) NLO + (F) pow (a s , a ). So 
if the next-to-leading (NLO) prediction {F)mlo 
exists, one can fit for as and for the universal 
parameter oo. 

In addition to the power corrections, next-to- 
leading log (NLL) resummations can help to im- 
prove the description of differential event shape 
distributions. The use of NLL sums requires a 
matching with the NLO prediction in order to 
avoid double-counting of terms. 

Fitting the means of the event shape variables 
with NLO+PC predictions one obtains the results 
for as and a~o shown in Figure ^ The resulting 
as values are consistent with the current world 
average, as — 0.118. The data suggest a® = 0.5 
consistently for all event shape variables except 
for the C parameter, the fit of which is extremely 
sensitive to the fit range in the ZEUS analysis (in 
contrast to HI measurements) . This result for oo 
is roughly compatible with the results from other 
experiments. 

3. PROMPT PHOTON PRODUCTION 

The ZEUS collaboration has for the first 
time measured prompt photon production in 
DIS in the full HERA 1 data sample of about 
120 pb -1 [2]. Prompt photons are required to be 
well-isolated and are identified using calorimeter 
cluster shape information; neutral particle back- 
grounds are subtracted using Monte Carlo (MC) 
predictions. 

Two analyses were applied to the data: In the 
first, inclusive prompt photon production is stud- 
ied, and the data are compared to the leading- 
order (LO) MC models HERWIG and PYTHIA; 
none of the models is able to describe all features 
of the data, and both are far away from the data 
in terms of normalisation. 

In the second analysis, a jet in the central 
rapidity region was required in addition to the 
isolated photon. Jets were reconstructed using 
a cone algorithm in the laboratory frame; their 
transverse energy had to be in excess of 6 GeV, 
and their pseudo-rapidity between -1.5 and 1.8. 
The measurements, which yield a total cross- 
section of 0.86 pb for prompt photon plus jets, 
are again compared to HERWIG and PYTHIA, 
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Figure 2. ZEUS prompt photon plus jet sam- 
ple: photon pseudo-rapidity distribution. The 
data, which are shown together with the energy 
scale uncertainty band, are compared to the NLO 
calculations for which the renormalisation scale 
uncertainty is also indicated. Also shown is the 
contribution to the cross-section which is due to 
prompt photon emission from the quark line. 



but there exists also an NLO QCD calculation 
for this signature. PYTHIA and HERWIG man- 
age to describe the shapes of the photon and jet 
transverse energy distributions, but they fail in 
describing both pseudo-rapidity distributions. In 
addition, their normalisation is again of by factors 
of about 2 (4) for PYTHIA (HERWIG). 

The NLO QCD calculation for the prompt pho- 
ton plus jet analysis predicts a cross-section of 
about 1.33 pb, which is supposed to be lowered 
by hadronisation corrections by about 30-40 %. 
The normalisation of the calculation is therefore 
compatible with the data. As an example, Fig- 
ure shows the prompt photon pseudo-rapidity 
distribution. The data are in agreement with the 
NLO predictions; however, the statistical uncer- 
tainty of the measurement is very large. The 
same is true for the jet transverse energy distri- 
bution. However, the calculations predict a too 
large cross-section at low photon transverse ener- 
gies and at high (forward) jet pseudo-rapidities. 
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Figure 3. ZEUS as measurements based on the 
mean integrated jet shape (ip(r — 0.5)} in DIS 
events as a function of the jet transverse energy. 
The data are shown together with the combined 
as value and the world average. 



4. NEWS FROM JET PHYSICS 

4.1. Determinations of a s from jet shapes 

The ZEUS collaboration has measured the sub- 
structure dependence of jet cross-sections in pho- 
toproduction and DIS with the aim of extract- 
ing the strong coupling parameter as EH- ZEUS 
used mainly the mean integrated jet shape (ip{r)} 
which is defined as: 
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Here, Nj ets is the total number of jets in the sam- 
ple, Ex(r) is the transverse energy within a cone 
of radius r around the jet axis, and El^ 1 is the 
transverse energy of the jet. The jet shape was 
measured, for example, as a function of the jet 
transverse energy in DIS events. For each data 
point, an as value is determined (Figure 0J; the 
resulting values are evolved to the mass of the 
Z, and the different values are combined into one 
a s (M z ) value: a s = 0.1176±g;gg^, where only 
the dominating theoretical uncertainty is given. 

4.2. Summary on as 

Figure 0] gives an overview of various as mea- 
surements performed by the HERA collabora- 
tions in jet events. The measurements are in per- 
fect agreement with the world average; the al- 
ready very good precision of the measurements, 
which has been increased over the past few years, 



will be even more increased once the HERA 2 
data are collected and analysed. 
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Figure 4. A summary of as measurements from 
HI and ZEUS jet analyses. 



4.3. Jets in global QCD fits 

Jet data from the HERA experiments can be 
used for more than just determinations of the 
strong coupling. Recently, the ZEUS collabora- 
tion has started to include jet data into their 
global QCD fits 0|, thus supplementing the usual 
i*2 data. The aim of using jet data is to im- 
prove the precision of the PDF measurements, 
especially the gluon, at high values of x. 

The use of jet data in the fits, however, requires 
a fast evaluation of the NLO predictions for the 
jet cross-sections. During the minimisation pro- 
cedure, this prediction has to be performed many 
thousand times. The usual typical 8 hours for 
sufficient statistics clearly is not acceptable here. 
The way around this problem is based on the 
following method. The jet cross-section can be 
written as the convolution of the hard partonic 
cross-section a with the PDF / Q : 

a = E Q E^° =1 a5 / dr)f a (r), (j, F ) ■ a n (r), fi F ) (2) 
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Figure 5. ZEUS global QCD fit. The gluon den- 
sity as a function of x is shown for various Q 2 
bins. Presented are the results of the ZEUS fit 
including jet data. 



Here, the first sum runs over all parton flavours. 
Defining a sufficiently finely binned grid in (77, /j,p) 
and assuming that the PDF / a (?7, (j.f) is flat in ev- 
ery bin separately, one can write the cross-section 
in one bin of the (j],/j,f) grid as: 

c(»7,Mf) ~ S„/ u ()j,mf)-^ Qs J dija n (77, (i F ) (3) 

where the integration is now only over the corre- 
sponding bin in 77. The summation of all pieces 
a(rj,HF) leads again to the total cross-section. 
The benefit of this method is that the coefficients 
(the integrals in Equation[3Jl have to be evaluated 
only once since they are independent of the PDFs; 



they can thus be stored in a grid and be convo- 
luted with any PDF provided from the outside, 
for example from the fitting program. The con- 
volution then takes only parts of a second, which 
allows for the use of this cross-section evaluation 
in fitting programs. 

The effect of the use of jet data can be viewed 
in Figure El which shows the results for the gluon 
density as achieved by the ZEUS fit using inclu- 
sive (F2) and jet data. The data are shown as 
function of x in different bins of Q 2 . A clear de- 
crease of the size of the uncertainty for the low-Q 2 
region is found; the effect is less striking for higher 
Q 2 values. 

For the future, ZEUS aims at including fur- 
ther jet data from photoproduction, heavy flavour 
production and also, if possible, from proton- 
antiproton collisions into the QCD fit. 

4.4. FROM PHOTOPRODUCTION TO 
DIS 

The transition region from photoproduction to 
DIS has since long been problematic in the the- 
oretical description. The HERA collaborations 
have published new results on this topic |5lfil7| . 

HI has analysed dijet cross-sections in the kine- 
matic region 5 < Q 2 < 100 GeV 2 0. The data, 
taken in the years 1996 and 1997, are compared to 
direct NLO calculations. Overall good agreement 
between data and predictions is observed. How- 
ever, comparing the data to NLO predictions for 
dijet events with small azimuthal separation be- 
tween the two hardest jets shows a clear excess 
of the data over the prediction. This discrepancy 
can be cured partly using a NLO three-jet cal- 
culation; however discrepancies remain at low Q 2 
and x. 

Figure shows the ZEUS dijet cross-section as 
a function of Q 2 over a wide kinematic range from 
GeV 2 to 2000 GeV 2 6 . The data are sepa- 
rated into regions of direct- and resolved-enriched 
events by means of the observable x 7 which quan- 
tifies the fraction of the photon's energy that 
entered the hard scattering. Shown is also the 
NLO QCD prediction. For the DIS regime, the 
NLO calculation, which does not include resolved 
contributions, clearly fails to describe the low-x 7 
component (independent of Q 2 ) and also the low- 
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Figure 6. ZEUS dijet cross-section as a func- 
tion of Q 2 . Separately shown are the direct- and 
resolved-enriched samples and their description 
by the NLO QCD calculation. 
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Figure 7. HI triple-differential dijet cross-section. 



Q 2 region. However, the prediction is very sensi- 
tive to the renormalisation scale choice, and scale 
choices can be made such that the data are com- 
patible with the predictions. For the photopro- 
duction region, which also is dominated by low 
Xj values, the (direct +resolved) NLO calculation 
is able to describe the data. 

HI has measured the triple-differential dijet 
cross-section d 3 a/dErdQ 2 dx 1 jj] {Et is the 
mean transverse energy of the two jets). Also 
here the conclusion is that at low Q 2 and also 
at low values of Et NLO QCD calculations fail 
to describe the data at low values of x 1 . This is 
the case even if the calculations include the con- 
tributions from the resolved photon (Figure 0) . 
In the same analysis, HI could also demonstrate 
the importance of the initial and final state par- 
ton showers which help considerably in describing 
the data with LO MC models. 

4.5. PARTON DYNAMICS AT LOW x 

Another problematic issue in jet physics is the 
question of parton dynamics at low x. Earlier 
studies of so-called forward jets |H| were designed 



to hunt for the break-down of DGLAP evolution 
and the onset of BFKL signatures. However, 
the results were never really conclusive, mainly 
because no NNLO DGLAP calculation and no 
BFKL prediction usable for experimentalists were 
(and are) available. Also the influence of resolved 
contributions to the predictions was unclear. 

There are new results by the HI collabora- 
tion on forward jet and ir° production in low 
Q 2 DIS prO], The forward jet analysis, which 
was performed in the kinematic region 5 < 
Q 2 < 85 GeV 2 , studied the triple-differential 
cross-section d 3 a/dp T dQ 2 dx and compared the 
data to NLO QCD calculations and various MC 
models [5] . One result is shown in Figure|Sl where 
the data are presented as a function of x and are 
compared to the NLO QCD calculation which 
clearly fails to describe the data at low values 
of x. MC models implementing resolved photon 
contributions and models using the the colour- 
dipole model come close to the data, although 
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Figure 8. HI forward jet measurement as a func- 
tion of x. The data are compared to (direct) NLO 
QCD calculations. Also indicated are the energy 
scale and renormalisation scale uncertainties. 



a discrepancy between data and predictions re- 
mains at lowest values of x. A more differential 
study shows that the discrepancies are largest for 
low transverse energies and low Q 2 . 

For the case of forward ir° production JU| , most 
of the available predictions do a sufficiently good 
job. Figure shows the forward ir° cross-section 
as a function of x in three bins of Q 2 ; the good 
description of the data by most models (including 
the NLO calculation involving a convolution of 
matrix elements and fragmentation functions) is 
prominent. 

5. STRANGE PENTAQUARKS 

The measurement of narrow baryonic reso- 
nances of positive strangeness in the K + n and 
Kgp decay channels with a mass of about 
1530 MeV and a very narrow width has recently 
raised a lot of attention. The findings are consis- 
tent with five-quark bound states, so-called pen- 
taquarks, + = uudds. 

The ZEUS collaboration started a search for 
the strange (anti)pentaquark in the K s p(p) chan- 
nel relying mainly on the precise tracking 
information from the central drift chamber. The 
analysis used all events with Q 2 > 1 GeV 2 of the 
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Figure 9. HI forward pion measurement as a 
function of x. The data are compared to vari- 
ous QCD calculations and models. 



full 1996-2000 data statistics of 121 pb -1 , about 
1.600.000 events. 

The Kgp(p) invariant mass was obtained by 
combining the K s and (anti)proton candidates 
and fixing the kaon mass to the PDG value. Fig- 
ure El shows the resulting pentaquark candidate 
mass distribution together with the fit (a back- 
ground function and two gaussians), the back- 
ground function and the prediction of the ARI- 
ADNE Monte Carlo. A clear peak at 1522 MeV 
with a width of about 6 MeV is visible. The num- 
ber of events ascribed to the signal is 221±48, of 
which 96±34 are in the antiproton channel. Thus, 
not only could the pentaquark be confirmed, but 
also the first observation of the antipentaquark 
with quark content uudds was made. 

The charmed pentaquark state that was ob- 
served by the HI collaboration will be discussed 
in another proceedings contribution |12|. 
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Figure 10. The ZEUS pK° s (pK° s ) invariant mass 
distribution with the peak at about 1522 MeV. 
Also shown are the results of a fit with a back- 
ground function and two gaussians, the back- 
ground function and the simulation of the ARI- 
ADNE Monte Carlo. 



between standard NC DIS events and instanton 
events using multi-dimensional cuts. Instanton 
events were simulated using the QCDINS calcu- 
lation which predicts an instanton cross-section 
of 8.9 pb. Under the assumption that all se- 
lected data events are instanton events, an upper 
limit on the instanton production cross-section of 
g < 26 pb was derived at 95 % C.L. 
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6. OTHER SELECTED TOPICS 

6.1. QCD instantons 

Instantons are non-perturbative fluctuations of 
the gauge fields in non-abelian gauge theories. 
But although instantons are predicted by the 
standard model, they have not been observed so 
far. Especially QCD instantons, however, are 
supposed to have short-distance implications at 
sufficiently low energies; they are, for example, 
supposed to induce characteristic topologies in 
DIS ep scattering events. Results on instanton 
searches have been reported on by the HI collab- 
oration before ^3]. Recently the ZEUS collab- 
oration has published the results of a search for 
instanton events |14| performed in about 40 pb -1 
of data collected in the years 1996/97. The data 
were selected according to Q 2 > 120 GeV 2 . For 
the search for instanton events, several discrimi- 
nating techniques were applied which distinguish 



Figure 11. Bose-Einstein correlations: ZEUS 
measurement of R(Q\2) for Q 2 > 4 GeV 2 to- 
gether with two fits. The data points included 
in the fit are marked with circles. 



6.2. Bose-Einstein correlations 

ZEUS has measured Bose-Einstein (BE) corre- 
lations in NC DIS over a wide kinematic range 
in Q 2 from 0.1 GeV 2 to 8000 GeV 2 PS] in over 
120 pb- 1 of data collected from 1996 to 2000. 
As a measure for the correlation, the analysis 
used the double ratio R(Q 12 ) = C, data /C, MC ' noBE ', 
where Q\2 is the Lorentz-invariant momentum 
difference of the particles (which are assumed to 
be pions), and ( data (C, MC ' noBE ) are the ratios of 
inclusive two-particle densities for like-sign and 

unlike-sign bosons, C = /?(++; )// 5 (^ — )i f° r 

data and Monte Carlo, respectively. The Monte 
Carlo sample was generated without simulating 
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Bose-Einstein effects, assuming that BE correla- 
tions can be factorised from other types of cor- 
relations and that non-BE correlations are well 
described by the MC models. Therefore R(Qi2) 
should be sensitive to BC correlations only. 

Figure ITT1 shows the double ratio R(Qi2) to- 
gether with two fit models which aim at extract- 
ing the size r of the source of final state pions 
and the coherence parameter A. The data points 
marked with circles were actually included in the 
fit; as expected the parameters extracted from 
the gaussian fit show no dependence on Q 2 . The 
radius of the production volume is found to be 
about r ~ 0.65 fm, and the coherence parame- 
ter A ~ 0.47. A two-dimensional study of the 
Bose-Einstein effect shows that the pion-emitting 
region is elongated. All findings are compatible 
with measurements from other experiments and 
colliders. 

6.3. Anti-Deuteron production 

The HI collaboration has measured anti- 
deuteron production |16| and compared the pro- 
duction rate to results from central proton-proton 
(Au-Au) collisions measured at the CERN ISR 
(RHIC). 

No significant signal can be found for positively 
charged tracks, but for negatively charged parti- 
cles a clear signal of 45 anti-deuterons with an es- 
timated background of about 1 event is observed, 
leading to a total anti-deuteron cross-section of 
about 2.7 nb. 

The production of nuclei in particle collisions 
can be described in terms of the coalescence 
model in which the cross-section a a for the pro- 
duction of a nucleus with A nucleons can be re- 
lated to the cross-section for the production of 
free nucleons in the same reaction, cr/v, by a so- 
called coalescence parameter Ba- 

The (anti-)deuteron coalescence parameter 
value B-2 ~ 0.01 found by ZEUS is compatible 
with values deduced at lower center-of-mass en- 
ergies in pp and pA collisions, but much larger 
than observed in Au — Au collisions at similar 
nucleon-nucleon center-of-mass energies. This is 
understandable if the size of the hard reaction 
'fireball' is much smaller in pp / jp than in heavy 
ion collisions. 



7. CONCLUSION 

Jet and hadronic final state physics at HERA 
has been reviewed. Clear progress especially in jet 
physics in terms of statistical precision has been 
made over the past few years. Jets have even been 
used in global QCD fits with the aim of improving 
the knowledge of the gluon density at high values 
of x. In addition, exotic signatures like instan- 
tons or pentaquarks have been analysed. Never- 
theless, HERA physicists are eagerly waiting for 
HERA 2 data which promise enough statistics to 
solve remaining questions and to come to firmer 
conclusions concerning as and the PDFs. 
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